The axial N to A transition form factors from Lattice QCD 



O 

o 

(N 



On 
(N 

> 
O 
m 
o 
l> 
o 
\o 
o 



X 



C. Alexandrou Th. Leontiou °, J. W. Negele \ and A. Tsapalis " 
" Department of Physics, University of Cyprus, CY-1678 Nicosia, Cyprus 
Center for Theoretical Physics, Laboratory for Nuclear Science and Department of Physics, 

Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, U.S.A. 
Institute for Accelerating Systems and Applications, University of Athens, Athens, Creece. 

(Dated: February 1, 2008) 

We evaluate the A'^ to A axial transition form factors in lattice QCD with no dynamical sea quarks, 
with two degenerate flavors of dynamical Wilson quarks, and using domain wall valence fermions 
with three flavors of staggered sea quarks. We predict the ratio C^{q^)/CY (q^) relevant for parity 
violating asymmetry experiments and verify the off-diagonal Goldberger-Treiman relation. 
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The electromagnetic structure of the nucleon including 
the transition form factors for electroproduction of the A 
has been the subject of recent experimental |l| and the- 
oretical studies [3, [1] ■ The N to A transition has the ad- 
vantages that the A(1232) is the dominant, clearly acces- 
sible nucleon resonance, and the isovector spin-flip tran- 
sition provides selective information on hadron structure. 
The weak structure functions provide valuable input of- 
ten complimentary to that obtained from electromagnetic 
probes. In the nucleon, for example, measurements of 
the elastic parity violating asymmetry yield information 
on strange quark contributions. The N to A transition 
filters out the isoscalar ss contributions, so the parity vi- 
olating asymmetry in weak neutral and charge changing 
N to A transitions probes isovector structure [1, Q not 
accessible in the study of strange isoscalar quark cur- 
rents. Furthermore its isovector nature could be used 
to detect physics beyond the standard model since it is 
sensitive to additional heavy particles not appearing in 
the standard model. However, in order to compete with 
other low-energy semi-leptonic measurements, the N to 
A transition requires a determination of the parity vio- 
lating asymmetry to significantly better than 1% preci- 
sion Hence, a lattice prediction for the axial form fac- 
tors provides valuable input for ongoing experiments Q. 
In this work, we evaluate the dominant contribution to 
the parity violating asymmetry, determined by the ratio 
C^/Cg'. This is the off-diagonal analogue of the gA/gv 
ratio extracted from neutron /3-decay and therefore tests 
low-energy consequences of chiral symmetry, such as the 
off-diagonal Goldberger-Treiman relation. In addition, 
the ratio of axial form factors C^/C^ provides a mea- 
sure of axial current conservation. 



Since this is the first lattice computation of the ax- 
ial N to A transition form factors, the starting point is 
an evaluation in the quenched theory using the standard 
Wilson action. A quenched calculation allows us to use 
a large lattice in order to minimize finite volume effects 
and obtain accurate results at small momentum trans- 
fers reaching pion mass, m^, down to about 410 MeV. 
In order to study the role of the pion cloud, which is 
expected to provide an important ingredient in the de- 
scription of the properties of the nucleon system, one 



requires dynamical configurations with light quarks on 
large volumes. In this work, the light quark regime is 
studied in two ways. First, we use configurations with the 
lightest available dynamical two fiavor Wilson fermions 
spanning approximately the same pion mass range as in 
the quenched calculation 0] - Second, we use a hybrid 
combination of domain wall valence quarks, which have 
chiral symmetry on the lattice, and MILC configurations 
generated with three flavors of staggered sea quarks us- 
ing the Asqtad improved action [8]. The effectiveness of 
this hybrid combination has recently been demonstrated 
in the successful precision calculation of the axial charge, 
gA 0- Since Wilson fermions have discretization errors 
in the lattice spacing, a, of 0{a) and break chiral sym- 
metry whereas the hybrid action has discretization er- 
rors of C(a^) and chirally symmetric valence fermions, 
agreement between calculations using these two lattice 
actions provides a non-trivial check of consistency of the 
lattice results. The hybrid calculation is the most com- 
putationally demanding, since it requires propagators on 
a five-dimensional lattice. The bare quark mass for the 
domain wall fermions, the size of the fifth dimension and 
the renormalization factors, Za, for the four-dimensional 
axial current are taken from Ref. [9|. As in the case of 
Wilson fermions, we consider three values of light quark 
mass with the strange sea quark mass fixed to approxi- 
mately its physical value Q. In all cases we use Wup- 
pertal smeared ^ interpolating fields at the source and 
sink. In the unquenched Wilson case, to minimize fiuc- 
tuations TT] we use hypercubic (HYP) smearing [T2| on 
the spatial links entering in the Wuppertal smearing func- 
tion at the source and sink whereas for the hybrid case 
all gauge links in the fermion action are HYP smeared. 
We list the parameters used in our computations in Ta- 
ble HI The value of the lattice spacing is determined from 
the nucleon mass at the physical limit for the case of 
Wilson fermions and for the staggered sea quark configu- 
rations, we take the value determined from heavy quark 
spectroscopy (l3j . 



The invariant N to A weak matrix element, expressed 
in terms of four transition form factors [l^, [3| can be 
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where = p' ^ ~Pfi is the momentum transfer, A^{x) = 
'ip{x)j^j5^ip{x) is the isovector part of the axial cur- 
rent, and is the third Pauh matrix. We evalu- 
ate this matrix element on the lattice by computing 
the nucleon two-point function {t;p;T4), the A two- 
point function Gf^(t; p; r4), and the three point function 
(G^^''^(t2,ii;p',p;r)) and forming the ratio [3 
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where the indices i are summed, r4 = i 
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'^^ Q J . For large time separations between ti , 

the time when a photon interacts with a quark, and t2, 
the time when the A is annihilated, the ratio of Eq. ^ 
becomes time independent and yields the transition ma- 
trix element of Eq. ([T]). The source-sink time separation 
is optimized as in Ref. [ll| so that a plateau is clearly 
identified when varying ti . We use kinematics where the 
A is produced at rest and = —q^ is the Euclidean 
momentum transfer squared. There are various choices 
for the Rarita-Schwinger spinor index, cr, and projec- 
tion matrices, F, that can be used in the computation 
of the three point function, each requiring a sequential 
inversion. We use this freedom to construct optimized A 
sources that maximize the number of lattice momentum 
vectors contributing to a given value of in analogy 
with the evaluation of the electromagnetic N to A tran- 
sition form factors 

In Fig. [T] we show quenched and unquenched results 
obtained with Wilson fermions. All errors shown are ob- 
tained using a jackknife analysis. We observe that 
is consistent with zero and that unquenching effects are 
small for the dominant form factors, and C^. In 
contrast, for the form factor C^, dynamical fermions 
produce a dramatic increase at low momentum trans- 
fer relative to the quenched results. Such large devia- 
tions between quenched and full QCD results for these 
relatively heavy quark masses are unusual, making this 
an interesting quantity with which to study effects of 
unquenching. In Fig. ^ we compare the Wilson and 
hybrid results for the two dominant form factors for 
ttItt ~ 500 MeV. The hybrid results thus corroborate 



TABLE I: The number of configurations, the hopping pa- 
rameter, K, for Wilson fermions or the light-quark mass, mi, 
for staggered quarks, and the pion, nucleon and A masses in 
lattice units. 



no. confs 


K or ami 


am,r 




oMa 


Quenched 32=* x 64 a''^ = 2.14(6) GeV 


200 
200 
200 


0.1554 
0.1558 
0.1562 
=0.1571 


0.263(2) 
0.229(2) 
0.192(2) 
0. 


0.592(5) 
0.556(6) 
0.518(6) 
0.439(4) 


0.687(7) 
0.666(8) 
0.646(9) 
0.598(6) 


Unquenched Wilson 24^ x 40 [6] a'^ = 2.56(10) GeV 


185 
157 


0.1575 
0.1580 


0.270(3) 
0.199(3) 


0.580(7) 
0.500(10) 


0.645(5) 
0.581(14) 


Unquenched Wilson 24=^ x 32 [7] a'^ = 2.56(10) GeV 


200 


0.15825 
K.C = 0.1585 


0.150(3) fl4] 
0. 


0.423(7) 
0.366(13) 


0.533(8) 
0.486(14) 


MILC 20^ X 64 a"^ = 1.58 Ge> 


/ 


150 
150 


0.03 
0.02 


0.373(3) 
0.306(3) 


0.886(7) 
0.800(10) 


1.057(14) 
0.992(16) 


MILC 28^ X 64 a"! = 1.58 Ge^ 


/ 


118 


0.01 


0.230(3) 


0.751(7) 


0.988(26) 



the small unquenching effects, and similar behavior is 
observed at the heavier and lighter quark masses. For 
Cf, also shown in Fig. [2l we observe the same large un- 
quenching effects observed for dynamical Wilsons. In ad- 
dition, although they are not shown, hybrid calculations 
yield C^- - 0. A dipole Ansatz C5^(0)/(l + Q'^/Mlf de- 
scribes the Q^-dependence of well, as shown by the 
curves in Fig. [51 yielding, at this pion mass, an axial mass 
Ma ~ 1-8(1) GeV. In this range of pion masses, we ob- 
served a weak quark mass dependence for Ma that how- 
ever needs to be checked for lighter quarks before com- 
paring to the experimental result of 1.28±0.10 GeV [l7| . 

In the chiral limit, axial current conservation leads 
to the relation C^^Q"^) = MfjCiiQ^)/Q^. In Fig. [3 
we show the ratio {Q"^ /Mj^) Ci{Q'^)/cf{Q'^) for Wil- 
son and domain wall fermions at the lightest quark mass 
available in each case. The expected value in the chiral 
limit for this ratio is one. For finite quark mass, the axial 
current is not conserved and for Wilson fermions chiral 
symmetry is broken, so that deviations from one are ex- 
pected. We observe that this ratio differs from unity at 
low but approaches unity at higher values of Q^. For 
the quenched case, where we have accurate results, a lin- 
ear extrapolation in yields values that are consistent 

with unity for ~ 0.5 GeV^. It is reassuring that this 
chiral restoration is seen on the lattice even for Wilson 
fermions, confirming that the discrete theory is correctly 
representing continuum physics. For finite pion mass, 
deviations from unity are expected to be proportional to 
TO^ / (Q^ -|- mj) for chiral lattice fermions. As can be seen 
by the dashed line in the figure, this form describes well 
the Q^-dependence of the results obtained with domain 
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XN,=0 m^=0.56 GeV 
ON,=0 m_,=0.49 GeV 
*Np=0 m__=0.41 GeV 
*N,=2 m_^=0.69 GeV, Re£. 6 
• N,=2 m^=0.51 GeV, Ref. 6 
□ N =2 m =0.38 GeV, Ref. 1. 
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FIG. 1: The axial form factors Ci,Ct, Ci and Ci as a func- 
tion of . We show quenched lattice results at k = 0.1554 
(crosses), at k = 0.1558 (open circles) and at k = 0.1562 (as- 
terisks) and unquenched Wilson results at k = 0.1575 (filled 
triangles), k = 1580 (filled circles) and k = 0.15825 (open 
squares). We use = 0.8 



wall fermions. 

For finite mass pions, partial conservation of the ax- 
ial current (9pAJ^(a;) — /^mj7r°(a;)) and pole domi- 
nance lead to the ofF-diagonal Goldberger-Treiman re- 
lation c^(g2) 



fTrgnNAiQ )/'2Mn whcrc g^^NAiQ ) 
is determined from the matrix element of the pseu- 
doscalar density < A~^\'ip{x)^5^^p{x)\p > and the 
pion decay constant, is calculated from the two- 

point function < 0|^4(x)|7r >. To relate the lattice pion 
matrix element to its physical value, we need the pseu- 
doscalar renormalization constant, Zp^ and we use for 
quenched [l^ and dynamical Wilson fermions [13] the 
value Zp{^'^a^ ~ 1)=0.5. In Fig. 2] we show the ratio 
fng-RNA/ {'^MjyC^) for Wilson fermions, which is almost 
independent of Q^, indicating similar dependence for 
g-wNA and C^. The ratio becomes consistent with unity 

for ^ 1 GeV^, in agreement with the off-diagonal 
Goldberger-Treiman relation. The large deviations seen 
for the lightest quark mass for dynamical fermions point 
to lattice artifacts that become more dominant for small 
quark masses. 

The presently unmeasured ratio /C^ is an interest- 
ing prediction of our lattice calculation. The form factor 
can be obtained from the electromagnetic N to A 
transition. Using our lattice results for the dipole and 
electric quadrupole Sachs factors, Qmi and Qe2 [HI, we 
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FIG. 2: The dominant axial form factors (top) and 
(middle) and as a function of , for similar pion mass, 
for quenched at k = 0.1558 (open circles ) and dynamical 
Wilson fermions at k = 0.1580 (filled circles) and in the hybrid 
approach at ami = 0.01 (squares with a cross). The curves 
are fits to the dipole form C^(0)/(1 + Q'^/M\f. 

1.2 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 




* N =0 m =0.41 GeV 

P tr 

^ Np=0 ni^=0 (linear extrap.) 
n N =2 m =0.38 GeV, Ref. 7 



O MILC m =0.34 GeV 



_1_ 



_L_ 



_L_ 



0.5 



1.0 



1.5 



2.0 



Q2 (GeV2) 



FIG. 3; The ratio {Q^/M%) Ci/C^ is shown versus in 
the quenched theory at k = 0.1562 (asterisks) and in the 
physical limit (filled triangles), for dynamical Wilson fermions 
at At = 0.15825 (open squares) and in the hybrid approach at 
ami = 0.01 (open circles). The short dashed line denotes 
the soft pion fit described in the text. 



extract using the relation 



3 MA(MAr + MA) 
2{Mn + M^Y + Q' 



{Qmi — Ge2) ■ (3) 



From the ratio C^/C^ shown in Fig. [51 we observe 
that quenched and unquenched results at m^r ~ 500 MeV 
agree within errors, and that the quenched data have a 
sufhciently weak mass dependence that when extrapo- 
lated to the physical pion mass, the results nearly coin- 
cide with the 500 MeV data within errors. Hence, it is 
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low energy consequences of chiral symmetry Q , this ratio 
provides the basis for a physical estimate of the parity vi- 
olating asymmetry. Under the assumptions that ~ 
and is suppressed as compared to C^, both of which 
are supported by our lattice results, the parity violating 
asymmetry can be shown to be proportional to this ra- 
tio [4]. Thus, our lattice results show that this ratio and, 
to a first approximation, the parity violating asymmetry 
is non-zero at = and increases for ^ 1.5 GeV^. 

In summary, we have provided a first lattice calculation 
of the axial transition form factors, which are to be mea- 
sured at Jefferson Lab [5] . The first conclusion is that 
is consistent with zero whereas is small but shows un- 



FIG. 4: The ratio Uq^na/ {2MnC^) as a function of Q^. ^^^^^y ^-^-^ sensitivity to unquenching effects. The two 
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FIG. 5: The ratio C^/Cl' as a function of for quenched 
QCD at K = 0.1558 (open circles) and at the physical pion 
mass (filled triangles), for dynamical Wilson fermions a,t k — 
0.1580 (filled circles) and at the physical pion mass (crosses) 
and for the hybrid action at am; = 0.02 (squares with cross). 

reasonable to use the extrapolated quenched results as 
a first estimate of this ratio, which is the off-diagonal 
analogue of gA/gv in the nucleon. In addition to testing 



dominant form factors are and Cq , which are related 
in the chiral limit by axial current conservation. The 
ratio {Q^ /Mj^^ Cq /C^ is shown to approach unity as 
the quark mass decreases, as expected from chiral sym- 
metry. For any quark mass, the strong coupling, gjrNA, 
and the axial form factor, C^, show similar depen- 
dence, and the off diagonal Goldberger-Treiman relation 
is reproduced as the quark mass decreases. The ratio of 
C^/C^ , which provides a first approximation to the the 
parity violating asymmetry, is predicted to be non-zero 
at = with a two- fold increase when ~ 1.5 GeV. 
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